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Abstract 



It is shown that the single-step periodic signal (periodic telegraph signal) can 
not produce coherent stochastic resonance for diffusion on a segment with 
one absorbing and one reflecting end points while the multi-step periodic 
signal does. The general features of this process are exihibited. The resonant 
frequency is found to decrease and the mean first passage time at resonant 
frequency increases linearly, as we increase the length of the medium. The 
cycle variable is shown to be the proper argument to express the first passage 
time density function at resonance. A formula for first passage time density 
function at resonance is derived in terms of two universal functions, which 
clearly isolates its dependence on the length of the medium. 

PACS number(s):05.40.+j 
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I. INTRODUCTION 



After the pioneering achievement of separation of large DNA molecules in gel medium by 
the application of uniform and time-dependent periodic electric field the mechanism 
of cooperative interplay between random noise and a deterministic periodic signal attracts 
considerable interests. It has been found that with this technique, large molecules in the 
size range 2 to 400 kb exihibit size-dependent mobilities. Similar ideas have also arisen in 
other types of chromatographic processes 0. 

The first passage time (FPT) is a useful tool to investigate the diffusive transport prop- 
erty in a medium. The theory of first passage time has been worked out in great detail for 
both infinite medium and explicitly time-independent diffusive processes However, for 

explicitly time-dependent processes and in finite medium an analytic closed form expressions 
are not available. In this respect also this problem attracts much attentions to the scientific 
community. 

The first analysis of this phenomena has been done for a random walk on a lattice 
numerically, and for a diffusive process in a continuous medium with periodic signal of small 
amplitude perturbatively [|^ . Their results indicate that the oscillating field can create a form 
of coherent motion capable of reducing the first passage time by a significant amount. This 
enhancement of the mobility of a particle in a diffusive medium by the application of proper 
oscillating field is known in the literatures [§,0 as coherent stochastic resonance (CSR). 

In order to investigate the reason for this cooperative behavior of random noise and 
deterministic periodic signal this problem has been formulated in much simpler terms by 
approximating the sinusoidal periodic signal by the telegraph signal p[ and subsequently 
it has been shown @,10] that the telegraph signal can not produce CSR. Finally, when 
the sinusiodal signal has been approximated as a multi-step periodic signal [l^, the non- 
monotonic behavior of mean first passage time (MFPT) with respect to the characteristic 
frequency of the periodic signal is recovered explaining the reason of CSR in the case of two 
absorbing boundaries explicitly. The general characteristics of the moments of first passage 
time density function (FPTDF) in their calculations |10] for continuous medium and with 
arbitrary amplitide of the periodic signal are found to be in agreement with the numerical 
simulation of the random walk model on a lattice [0]. 

When the phenomenon of CSR is being discussed for the linear systems, it turns out that 
the boundary conditions play a crucial role in some models [p!I| -[I^. Doering and Gadoua 
have considered the jumps in a linear double- well potential when the potential fluctuates 
between two values at a rate 7. The non-monotonic dependence of the MFPT on 7 (called " 
the resonance activation" ) has been found to occur when one of the boundaries is absorbing 
and the other is reflecting Brey and Cassado-Pascual |T2I] have investigated random 



walk on a one dimensional lattice in which at any time each site has one of two transition 
rates which are being allowed to change at random times. Similar effects of boundary condi- 
tions have been observed. Linear models with asymmetric ||13[ and symmetric |jl^ random 
telegraph signals (dichotomous noise) also show that MFPT behaves non-monotonically with 
the jump rates for the transition between these two states when one of the boundaries is 
absorbing and the other is reflecting. 

In this paper we consider an overdamped linear system driven by a sinusoidal force 
and embedded in a noisy environment which is taken to be Gaussian white. As mentioned 
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before, no analytic closed form expression is available in the literature for this explicitly time 
dependent problem. We therefore approximate the sinusoidal force (signal) by a multi-step 
periodic signal (explained below). This is the system which we considered before |Ty]. The 
only difference is that in this paper we have asymmetric boundary conditions, i.e., one end 



point is a reflecting boundary and the other one is absorbing while in the previous work |T0 
both the boundaries were taken as absorbing. 

We note that the linear system with single-step periodic signal (periodic telegraph signal) 
could be recovered as a special case of this system with multi-step periodic signal. The linear 
system with single-step periodic signal (periodic telegraph signal) with one reflecting and one 
absorbing boundaries has been considered in the literature . Non-monotonic behavior of 
MFPT with respect to the characteristic frequency of the periodic telegraph signal has been 
reported In this paper we show that the telegraph signal does not produce any non- 



monotonic behavior of MFPT with respect to frequency. This result is in direct contradiction 
with the result obtained in . The explanation of this contradiction is stated in the text 
below (in Sec. II and Sec.IIIA). 

The paper is organised as follows. In Sec. II, we give the formulation of the problem. 
The basic structure of this formulation is similar as before ||10|]. The only change with the 
previous one is to incorporate the effect of asymmetric boundary conditions. For the sake 
of completeness we rewrite the final formulae here with incorporation of proper boundary 
conditions. After giving the derivation of MFPT in Sec. II, the results of the calculations 
are discussed in Sec. III. First we present the general characteristics of CSR. The calculation 
clearly exhibits how resonance appears in our multi-step approximation and fails to show in 
single-step telegraph approximation of the periodic signal. The general characteristics of the 
moments and the characteristic features of FPTDF for this phenomena are also presented 
in this subsection. In the next subsection we focus on the resonance point and demonstrate 
some special features associated with it. In particular, we show in this subsection (Sec.IIIB) 
that the cycle variable is the proper argument to express the FPTDF at resonance. Fur- 
ther, it is shown that the FPTDF at resonance can be expressed in terms of two universal 
functions. This feature clearly isolates its dependence on the length of the medium. Finally, 
few concluding remarks have been added in Sec. IV. 



II. FORMULATION OF THE PROBLEM 

We consider diffusion in one dimension perturbed by a periodic force. The motion of the 
particle is given by the Langevin equation 

X = Asmnt + ^{t), (1) 

where X refers to the stochastic variable, A and Q are the amplitude and frequency of the 
sinusoidal signal and ^(t) is a zero mean Gaussian white noise of strength D with auto- 
correlation function given by 

<mat')>=^mt-t')- (2) 

The motion is confined between a reflecting boundary at x = and an absorbing boundary 
at X = L. The Fokker Planck equation corresponding to Eq.(|l]) is 
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(3) 



where p{x,t) and j{x,t) refer to the probabihty density and probabihty current density 
respectively at position x and at time t. The reflecting boundary condition at x = implies 
that j{0,t) = and absorbing boundary conditions at x = L suggests that p{L,t) = 0. We 
now introduce the dimensionless variables 

e = {A/D)x, e = {A^/D)t, uj = Q/{AyD), (4) 

to write Eq.(^) in terms of new variables: 

The boundary conditions are rewritten as j{0,6) = p{A,6) = 0, where A = {A/D)L. In 
the following we calculate all the physical quantities in terms of these new variables and 
if required, one may translate all the interpretations in terms of the usual variables by the 
transformation equations Eq.^. 

No analytic solution exists for Eq. (P) with the boundary conditions mentioned. We thus 
introduce a scheme to approximate the force sinuO as a multi-step periodic signal []TU[. 



This scheme is in contrast to the procedure adopted in where they discretise the Eq.(^ 
using finite difference method and simulate the problem on a lattice of space and time. The 
construction of multi-step periodic signal is as follows. We divide the half cycle of the signal 
by {2p + 1) intervals so that each interval in the horizontal ^-axis is of size {A9/{2p + 1)) 
with uj A 9 = TT. We define {2p + l)numbers along the vertical axis as 

Sk = — — ■,k = l,2,...,p (6a) 

Sp+i = 1 (6b) 
Sp+i+r = Sp+i-r ; r = 1, 2, p. (6c) 

Each number is associated with the interval ''^2p+f'^ < ^ < with /c = 1, 2, (2p + 1). 
The Eq.(^ clearly shows that 

< Si < S2 < ... < Sp < Sp+i = 1 > Sp+2 > Sp+3 > ... > S2p+1 > . (7) 

Eq.(|3) states that in order to reach the maximum value (= 1) of the signal from the zero 
level we have to have {p + 1) step up and from the maximum to the zero level we have 
{p + 1) step down. This is for the positive half-cycle. For the negative half-cycle similar 
constructions have been done with the replacement Sk —Sk,\/k and each number —Sk is 
associated with the interval A^[l + ^] < 6 < Ae[l + ^] with = 1, 2, {2p + 1). 
This aproximation for the full one cycle of the sinusoidal signal (as shown in Fig.jlDis then 
repeated for the next successive cycles. The construction clearly shows that we get back 
the usual telegraph signal with p = 0. Approximation of the sinusoidal signal by the usual 
telegraph signal has been made by in the case of two absorbing boundaries. 

One may, however, note that the uj which we have defined for this approximated signal 
is not the same as that of sinusoidal signal, because the Fourier transform of sinusoidal 
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signal would give only one frequency while this approximated signal in the Fourier space 
corresponds to many sinusoidal frequencies, especially because of its sharp discontinuities. 
Yet we urge this approximation because in each small interval the equation becomes time- 
independent. 

The Fokker-Planck equation (Eq.(^) in each interval with this scheme will be that for a 
constant bias, namely 

where U\^) = —s{s > 0) for the positive half cycle and equals to +s{s > 0) for the negative 
half cycle with s being the value of Sk (Eq.(|^)) for the corresponding time interval. We wish 
to express the conditional probability density in each interval in terms of complete set of 
normalised eigenfunctions satisfying the boundary conditions mentioned above. For that it 
is convenient to cast Eq.(|D into an eigenvalue problem of Schrodinger type by setting 

9) = e-^'e^'^/'m . (9) 

Substituting the ansatz Eq.(|) in Eq.(|^) we obtain 

A = /x^ + sV4 , (10) 



0. 



(11) 



Reflecting boundary condition at ^ = 0, j{0,9) = 0, and absorbing boundary condition at 
^ = A, p(A, 9) = 0, associated with Eq.(H) with the help of Eq.(y) take the following form 



[-t/'0 + 0'k=o = O,0(A) = O. 



(12) 



In the future development we associate the index n for the positive half-cycle and index 
m for the negative. Index i will refer to the cycle number. 

Employing these boundary conditions (Eq.(^)) in some interval in the positive half- 
cycle, the normalised eigenfunctions, 0n(O ^^^e obtained as 



^ sin2/JnA 

2/ln 



1/2 



sin/i„(A - ^) 



(13) 



where obtained by solving the transcendental equation 

sin /i„A H — — cos /i^A = , 



(14) 



with s being the value of Sk (Eq.(P)) for the corresponding time interval where the conditional 
probability density is to be decomposed. The corresponding eigenvalues A„ are given by 



(15) 
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We note that A in Eq.@ must be positive. Therefore in the range — < /i^ < 0, 
employing the above procedure one obtains the eigenfunction 0(^) ~ sinh — A), where k 
would be obtained as a solution of the transcendental equation 



sA sinh kA + 2kA cosh kA = 



(16) 



where h? = —jJ^^K > 0. As no solution exists for Eq.(|T6D, we say that for the positive 
half-cycle the normalised eigenfunctions (pniO Eq. ([TB|) form a complete set. On the other 
hand, for the negative half-cycle, employing the similar procedure we find the normalised 
eigenfunctions as 



0m (0 



A 



sin 2/im A 



1/2 



sin/i„(A - 



(17) 



where /i^ are obtained by solving the transcendental equation 

sin /imA cos /imA = , 



(18) 



with s being the value of Sk for the corresponding time interval where the conditional prob- 
ability density is to be decomposed. We note that Eq.(p!7|) and Eq.(|18|) are true for /i^ 7^ 0. 
For = 0, boundary conditions suggest that one nontrivial eigenfunction exists only when 
sA = 2. Thus if the value of A is such that sA = 2, apart from the eigenfunctions given by 
Eq.(|l7D, one more eigenfunction exists which is 



^Mm=o(0 ( ^3 J 



A). 



(19) 



Further, in the range — < < 0, only one nontrivial eigenfunction exists for sA > 2. 
This is given by 



MO 



sinh 2kA 
2k 



-A 



1/2 



sinh /t(A — C.) 



(20) 



where = — /i^, n > and it is obtained by solving the transcendental equation 



— sinh K,A H cosh kA 

s 



0, 



while in the range < sA < 2, no nontrivial eigenfunction exists, 
eigenvalues in all the cases are given by 



2 



+ s'/A. 



(21) 

The corresponding 
(22) 



We thus see that the set of eigenvalues for a given s are different in positive and negative 
half-cycles because the transcendental equations, Eq.(|^) and Eq. ([T8|) , are different. Thus 
the corresponding eigenfunctions of Eq.(p!3[) and Eq.(|T^) are different functions, although 
they have the same form. Apart from these differences, we have seen that in the range 
— < /i^ < 0, there is always a nontrivial eigenvalue obtained from Eqs.(pT])-(^) (for 
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sA > 2) or (for sA = 2) and corresponding nontrivial eigenfunction given by Eq.(pOD 
( for sA > 2), or Eq. ([T9|) (for sA = 2 ) in the negative half-cycle while in the positive 
half-cycle, in this particular range, we do not have any nontrivial eigenfunction. For a given 
value of A, the corresponding set of {0n(O} forms the complete set of eigenfunctions in the 
positive half-cycle while the set {0m(O}(i^c^^'^i^S Eq.(|20|) or Eq.(|l^) as the case may be) 
forms the complete set in the negative half-cycle. 

As argued before, from our multi-step periodic signal the single-step periodic signal 
(telegraph signal) could be recovered with p = 0. In this situation s takes only one value, 
namely s = 1. The results given above for the eigenfunctions and eigenvalues are true 
for s = 1 (telegraph signal) also. As mentioned in the introduction this special case has 
been treated by Gitterman |T^. However the eigenfunctions and eigenvalues used there 
[see Eq.(14) in the paper |T^] are wrong because the eigenvalues do not satisfy the proper 
boundary conditions Eq.(0) for both the half-cycles. 

Once the complete set of normalised eigenfunctions and eigenvalues are determined the 
conditional probability density function p(^, 9 \ 0') in the positive half-cycle could be 
expressed as 

p{i.o\ i\e') = Y.<{i)u-{O^M-K{e~&)] , (23) 

n 

where 

«±(O = exp(±</2)0„(O . (24) 

with s as the corresponding value of Sk in the appropriate time interval where the condi- 
tional probability is being decomposed and 0„(O and A„ are given by Eq.(|I^) and Eq.(|T3p 
respectively. The conditional probability density function in any interval, say Z, can then be 
calculated from the previous history by convoluting it in each previous intervals: 




For the negative half-cycle the calculation of probability density function is similar except 
that we have to replace the index n by m and the probability density function is decomposed 

as 

v{i.e\i',o') = Y.u-^{i>UO^M-^m{e-e')] , (26) 

m 

where the expressions for ^^(0 same as in Eq . (p^ with n replaced by m and Am are 
given by Eq. (p^) . 

The survival probability at time 9 when the particle is known to start from ^ = at 
6^ = is defined as 

S{e\i,)= r«e,^Uo,0) . (27) 

JO 

The first passage time density function (FPTDF) g{6) is defined as 
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Physically, g[6)d6 gives the probability that the particle arrives at the absorbing boundary 
in the time interval 9 and 9 + d9. From this density function one can calculate various 
moments: 

roo 

< 9^ >= d99^g{9) . (29) 

From Eq.(pP|) one can easily calculate mean first passage time(MFPT) < 9 > and the 
variance o"^ =< 9"^ > — < 9 of the density function g{9). 

It is then quite straight-forward to calculate the survival probability at any interval of 
any cycle. We will write down the final formulae: 

^tki^ I ^o) = C'n(2p+i)(,_i)+fe 

xexp[-A„^,^^,^^^_,,^,(e-2(^-l)A0)] 

; [2(i - 1) + -^^1 A9<9< \2(i - 1) + —^1 A 9 , 

■,k = l,2,...,{2p+l) , (30a) 

SiA^ I ^o) = C'm(2p+i){,_i)+fe 

xexp[-A„^,^^,,^,_^,^,(^-(2z-l)A0)] 

^-^fe('"m(2p+i)(i_i)+fc) ' 

' [^2^ - ^) + ^] A < < [(2^ - 1) + ^] A , 

;A; = l,2,...,(2p+l) , (30b) 

where 

C:= rrfe«+(0 , (31a) 
Jo 

C-= l^d^u-^O , (31b) 
Jo 

and the functions Ok, are generated through the recursion relations: 

Oi(u- ) = Fi_i{u- ) , (32a) 

^'c('^n(2p+i)(i_i)+fc) ~ ^ '^n(2p+i)(i_i)+fc I ''^n(2p+i)(i_i)+(fc_i) ^ 

. (fc-1) Ag _ X ^1 

^^-^Pi 2p + 1 V'^"(2p+l){«-l)+fe '^n(2p+l){i-l) + (fe-l)JJ 
« J-V n(2p+i)(i_i)+(fe_i) / ' 

;A; = 2,3,...,(2p+l) , (32b) 
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-^^ ^(2p+l)(i-l) + l-' »T*(2p+l)(i-l) + l I "(2p+l)i 



xexp[- A0A„(,^^^,J 

x02p+i{u-^^^^J , (32c) 



r (A:-l) _ 

^ ^-^Pi 2p+l ^^™(2p+l)(i-i)+fc '^m(2p+i)(i_i)+(fc_i) jj 
^-^^-l('"m(2p+l){«-l) + (fc-l)) ' 

;A; = 2,3,...,(2p+l) , (32d) 



Fi(u„ ) = < u„ I ti™ > 

'V ri(2p+i)i+i/ n(2p+i)i+i I "i(2p+i)i 

X exp[- A eXm^^^^J 

x^2p+i(m+,,^^,)^) , (32e) 

with Fq{u~J = u~^{^q). The angular bracket in any equation imphes dot product of the 
corresponding functions, for e.g., 

<u+\u->= ^d^u+iO^^-iO ■ (33) 
Jo 

The cycle variable i runs over positive integers; i.e., z = 1, 2, 3, .... The positive and negative 
symbols of the survival probabilities indicate its value over positive and negative part of the 
cycles respectively. In all these expressions, viz., Eqs.(|30| ) - (|32|) , any subscript either n or m 
or both whereever they appear more than once the summation over them are implied. The 
effect of history is explicit in the expressions for survival probabilities. Once the survival 
probability S{6 \ ^q) is obtained from these formluae, the FPTDF, MFPT and the corre- 



sponding variance are obtained by employing Eqs.(p^)-([29|). Evaluation of MFPT and other 
relevant quantities requires sum of infinite series which must be truncated in order to obtain 
a final result. Convergence of MFPT is ensured by gradually increasing the number of terms 
(i.e., number of eigenvalues) for the calculation. The process is truncated when MFPT does 
not change upto two decimal point of accuracy with the change of number of terms. 

III. RESULTS AND DISCUSSIONS 

The survival probability, mean first passage time(MFPT), corresponding variances and 
first passage time density functions(FPTDF) are calculated using the derived formulae for 
this process. The results are summarised below. 

A. General features of CSR 

The MFPT is calculated for single-step telegraph signal {p = 0) with = A/2. Most 
of the calculations are done with this specific value of ^o- The variation of the results with 
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variation of is also demonstrated [see the text below]. No nonmonotonous behavior is 
observed in MFPT as we vary the frequency uj. This is in complete disagreement with 



Gitterman's observation |T^. Because of the wrong set of eigenvalues and corresponding 
eigenfunctions in both positive and negative half-cycles used to express the conditional prob- 
ability density, Gitterman observed nonmonotonic dependence of MFPT with respect 
to the characteristic frequency of the periodic telegraph signal. 

The calculation is done for the length A = 20 and the result is shown in the curve (a) 
of Fig. 2. However, when we take p = 1, i.e., when the sinusodal signal is approximated 
by two-step periodic signal, the calculation of MFPT for the same length shows clearly the 
nonmonotonous behavior. This is shown in curve (6) of the same figure. This result clearly 
demonstrates that mere flipping of the bias (signal) direction periodically would not produce 
the coherent motion. As the rate of flipping increases it merely prevents the particle more 
to reach the absorbing boundary and therefore MFPT increases monotonically. It may be 
noted that when the flipping rate is very high, the effect of signal is almost nil and the 
transport is effectively diffusive in nature. This is of course true in any type of periodic 
signal. Therefore, for any type of approximation of the sinusoidal signal or for any value of 
p, this feature would show up. In particular, for p = 1, we observe from curve (6) of fig. 2 
that MFPT would asymptotically reach the diffusive limit 3A^/8 (= 150 in this case). The 
usual telegraph signal offers a constant bias of maximum magnitude for the larger time than 
for a two-step approximation. Hence the particle always has a larger probability of reaching 
the absorbing boundary in short time for p = case than for p > case. Hence MFPT for 
p = and for any uo is always less than for p > case. This is observed in Fig. 2. 

The application of any bias always reduces the MFPT than for the non-biased diffusion. 
In GSR we always have a competition between diffusion and oscillatory effect of the bias. 
For very large frequency as the bias effect becomes ineffective MFPT would essentially be 
guided by diffusive process. For zero frequency of the multi-step periodic signal the MFPT 
can be analytically evaluated. When it starts from the mid-point of the medium it expresses 
as 

/2g-3siA/4\ 

< e{uo = 0, ^0 = A/2) >= 0.5(A/si) - ( j sinh(siA/4) . (34) 

When frequency is very small, the process is predominantly diffusion with constant value 
Si = 0.5sin(2^^) effective for < 6^ < ^2p+Ty- However, as frequency increases slowly, the 
probability of having increased bias value S2 (=1 for p = 1) before it reaches the absorbing 
boundary increases. This bias force reduces the survival probability and also MFPT. Hence 
one would expect a minimum to MFPT. On the otherhand, for usual telegraph signal (p = 
case), for very low frequency, from the very beginning bias force affects the particle with its 
maximum strength. When the frequency is very low, this constant bias diffusion continues 
for a longer time and there is no change-over of the magnitude of the bias as in the case of 
p = 1. After having a flip, the particle again suffers a constant bias diffusion in the direction 
opposite to the previous one and towards the reflecting boundary. As frequency increases 
slowly, this picture remains unchanged until a stage reaches for which the flipping effect 
becomes dominant during the particle's survivality inside the medium and MFPT increases. 
This is observed in Fig. 2. 

Next we continue all our calculation with p = 5 or, with six-step telegraph signal. Gal- 
culation reveals that the value of MFPT does not change much from that with p = 1. On 
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the other hand, p = 5 signal approximates the sinusidal signal better than p = 1 signal. We 
restrict our calculation with p = 5 approximation of the periodic signal. 

We calculate the MFPT < 6 > and the variance o"^ as a function of frequency u for 
different lengths (A = 10, 20, 30, 40, 50). These are presented in Fig. 3 and Fig. 4 respectively. 
Both the cumulants go through a minimum as frequency rises from very low value for each 
length A. It is observed that the minimum for both the moments occur almost at the 
same frequency for each length. This shows that the maximum cooperation between the 
deterministic signal and random noise occurs at this resonant frequency. The value of MFPT 
< 6 > increases with the length at all frequencies. This is understandable because as length 
increases, on an average the particle will spend more time in the medium before reaching the 
absorbing boundary. It is also observed that the frequency at which the minimum occurs 
shift towards low frequency as the length increases. 

Fig. 4 demonstrates the lowering of the dispersion at resonant frequencies confirming that 
the cooperation is maximum at these frequencies. Dispersion is more for higher lengths and 
as seen from the figure the dispersion merges to a specific value at very low frequency at 
various lengths. 

All the previous calculations are done when the particle starts initially from the mid 
point of the medium, i.e., .^o iii Eq.(^) is taken as A/2. At the length A = 30 the resonant 
frequency is found to be 0.06. The calculations are done one at resonant frequency and other 
two at the off- resonant frequencies {uj = 0.1 and uj = 0.0) when the particle starts from 
^0 = /5A where f3 lies between and 1. For zero frequency the MFPT can be analytically 
obtained. Its expression reads as 

A(l-8) /2e~^i^(i+^)/2\ 

< e{uj = 0, f3A) >= ^ - — . sinh(siA(l - (3)/2) . (35) 

si V si / 

The curves are shown in Fig. 5. It is evident that the value of < > is less for resonant 
frequency (curve (a)) than for its value for off-resonant frequencies(curves [b) and (c)). For 
each curve the maximum value of < > occurs at lower values of /3 or, when the particle 
starts from the left of the interval (near the reflecting boundary). Our signal starts with 
positive half-cycle and therefore the survival time of the particle would be more if the particle 
starts from the left of the interval. On the other hand, if it starts from right half of the 
medium (near the absorbing boundary), the initial surge of the signal helps the particle to 
reach the absorbing boundary more quickly. Hence average time of duration decreases. For 
curve (b) where the frequency of the signal {u = 0.1) is more than the resonant frequency, 
the oscillatory contribution is more than that for the resonant frequency making the MFPT 
large than for curve (a). For curve (c) where the frequency is zero {u = 0.0), the particle 
experiences a constant bias {s = si) towards the absorbing boundary all the time. This is 
in contrast to the other two cases where the particle has a probability to experience bias of 
magnitude more than si. As the value of si (the first value of six-step periodic signal) is 
smallest hence the curve (c) lies always above the curves (a) or (6). 

The FPTDF g{6] u) for different frequencies are calculated and plotted as a function of 
6 for a given length A = 30 [Figs. 6a & 6b]. One can clearly see the evolution of the FPTDF 
profile as the frequency increases from u = 0.001 to u = 0.1. 

At very small frequency (e.g. u = 0.001), the particle is acted on by a constant force 
s = Si (in the multistep periodic approximation) almost all the time before it reaches the 
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absorbing boundary. The FPTDF curve shows that we have only one maximum in the entire 
^-range. 

As the frequency slowly increases, one finds that, apart from only one maximum (like 
the one at very small frequency) other small peaks at larger 9 also show up gradually 
Also, these new small peaks start becoming stronger as the frequency increases. This is 
because the particle which initially sees a constant bias (s = si) for some time encounters 
an increased bias after a while (larger 6) and so on before reaching the absorbing boundary. 
As / g{9)d0 = 1, the area under the major profile decreases and is compensated by extra 
peaks at higher 9. We therefore see that as the frequency increases continuously, the small 
peaks that show up at higher 9 start growing and at the same time the first peak slowly 
decreases keeping the total area same. 

The resonance occurs at uj = 0.06 for A = 30, as shown in Fig. 6b. The MFPT has a 
minimum at that frequency. It has been argued that the synchronization between the signal 
and the noise is maximum at this frequency causing the enhancement of the probability of 
reaching the absorbing boundary at a short time. The FPTDF profile at resonance (at this 
frequency) evidently shows a dominant single peak. We also see that the constituent profiles 
of FPTDF adjust themselves in a very distinct way, as the frequency increases, in order to 
produce a dominant maximum at the resonance. 

Beyond the resonance (at higher frequencies) the peaks are quite numerous, distinct and 
identified separately, while the area of the dominant peak (at resonance) starts decreasing. 
In other words, the degree of synchronization is getting reduced as we go beyond the resonant 
frequency - if we identify the area of the peaks as the degree of synchronization (considering 
total area is normalized). 



B. Special features at resonance 

In this subsection we concentrate on the behaviour of the system at the resonance point. 
We have already discussed some general characteristics of CSR in the previous subsection. 
We find that for each length, A, a corresponding frequency uj* exists for which < 9 > and 
cr^ become minimum implying that the maximum cooperation between the deterministic 
periodic signal and random noise of the environment is taking place in helping the particle 
to reach the absorbing boundary. One therefore would naturally inquire about the relation 
of a;* with A. The curve of a;* as a function of A is plotted in Fig. 7. In the range of A we 
studied this curve is very well fitted with the formula 

cu* = C/A^ , (36) 

where C = 0.8053, and 7 = 0.7615. 

The values of MFPT at resonance < 9{u!*) > is plotted against the length A in Fig. 8 
and within the range of A we consider the relation between them is fitted to 

< 9{u*) >= 0.79A + 1.85 . (37) 

Of course, there will be deviation from this linear behaviour as A decreases further because 
< 9 > can not become positive for A = (corresponding to L = 0); < 9 > should be zero at 
A = 0. 
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We have already seen that at the resonance frequency we have single dominant peak of 
FPTDF, (7(^)[Fig.6b]. Since it is a general feature, for each length A we should get such 
behaviour. When we plot g{9)/uj* as a function of [a;*(A)^], we find that the curves for 
different lengths superpose over each other [Fig. 9] and the pattern of g{9)/uj* for different 
A or uj* is very similar, i.e., at particular values of [oJ*9], all curves show their maxima, and 
change in the behavioral patterns of the curves occur exactly at the same places of[co'*0]. 
This shows that [uj*0\ or the cycle number is the correct variable to describe the resonance 
behaviour. We may further note that such scaling of FPTDF would not be possible for 
any frequency other than the resonant frequencies because any frequency which is not the 
resonant frequency for one length may turn out to be the resonant frequency for some other 
length and the features of FPTDF are different for resonant and off-resonant frequencies as 
has been observed from Fig.Ga - Fig. 6b. The major dominant peaks of FPTDF g{9)/uj* for 
different lengths(A = 20,25,30,35,40,50) are drawn as a function of [oJ*9] in Fig. 9. The 
lowermost curve is for A = 20 and as length increases the upper curves are generated. The 
peaks for all the curves occur nearly at a quarter of a cycle. 

Having found the proper scaling of the argument of FPTDF, it is natural to enquire 
whether the FPTDF, / = g{9)/uj* can also be scaled properly, so that once f{x) with 
X = uj*9 is found for one length, the function f{x) can be obtained for any arbitrary length. 

In order to investigate this issue, first we observe that going from /(x;Ai) to /(a;;A2), 
where Ai and A2 are two different lengths we have to multiply one function by different 
amount depending on the value of the argument. Thus if at all any scaling relation exists, 
this should be of the form 

f,K{x) = /A(x)/i°(^') (38) 

It is indeed found to be true with a{x) given in Fig. 10. We note that a{x) is universal. 
From the relation (38) it is obvious that 

A-°(^)/a(x) = [M]-"^"Vm(^)- (39) 

The Eq.(^) shows that the function A~"'^^)/a(x) is independent of A and it is again universal. 
We verify this result in our calculation of FPTDF for different lengths and call this function 
U{x). The lnU{x) is also plotted in Fig. 10. Thus we could express /a (a;) at resonance in 
terms of two universal functions U{x) and a{x) as 

/A(a;) = f/(a;)A"(^) (40) 

The expression (40) clearly isolates the dependence of FPTDF on length of the medium A 
at resonance. The form (|40D for /a (a;) is plotted for a typical length A = 30 and compared 
with actual data points in Fig. 11. It is giving an excellent agreement. Further with this 
form, the MFPT is found for different A. It also yields required behaviour as in Eq. (p^) . 
This shows that not only we have found a scaling behaviour given by Eq.(^) which takes 
f{x) from one length to the other, it is also possible to obtain a simple expression Eq.(^OD 
for /a(x) for any arbitrary length. 
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IV. CONCLUDING REMARKS 



We consider a diffusive transport process perturbed by a periodic signai in continuous 
one dimensional medium fiaving one reflecting and one absorbing boundaries. We showed 
explicitly that the cooperative behaviour between the deterministic periodic signal and ran- 
dom noise leading to coherent motion occurs when the time- dependent sinusoidal signal is 
approximated by a multistep periodic signal and not with single-step telegraph signal. 

Although we study the process with six-step periodic signal, the formulation is quite gen- 
eral and applicable for any approximation with arbitrary number of steps. This formulation 
can also be applied to any arbitrary continuous periodic signal. Further, no perturbation of 
the signal amplitude is assumed in this formulation. 

It is observed that for large time oscillation of the signal plays a dominant role in the 
transport while in the low time regime frequency dependent bias force (i.e., the chance of 
having increased values of the bias in the same direction is more for high frequency than for 
low frequency) has the key factor. For very high frequency the bias effect is practically absent 
and the motion is purely diffusive in nature. At the resonance the maximum cooperation 
between the noise and the periodic signal takes place making the MFPT a minimum. 

The resonant frequency decreases and the mean first passage time at resonant frequency 
increases hnearly, as we increase the length of the medium. The corresponding relations are 
obtained. This fact provides a measure of the time scale of this process at resonance. 

Other important characteristic that we observe is that at the resonance the FPTDF for 
various lengths have similar behaviour as a function of cycle number. The curve shows that 
there is single dominant peak, which is a reflection of the fact that the synchronization 
between the deterministic periodic signal and the random noise is maximum at resonance. 
The peak positions of these curves occur very near to a quarter of a cycle showing that the 
cycle number is the correct argument to describe FPTDF at resonance. 

We also show that there exists a scaling relation between FPTDF at various lengths 
through some universal function a{x). The exact expression for FPTDF at resonance is 
obtained in terms of two universal functions, which clearly isolates its dependence on the 
length of the medium. This is a special feature of coherent stochastic resonance. This form 
may be of use in obtaining quicker result in cases where more complex situation is called 
for. 

There is also slight discrepancy in the position of the minimum of o"^ in comparison to 
the minima of < ^ > [Fig. 3, Fig. 4]. This may be due to the fact that all calculations are 
made to an end when the survival probability takes a value 1 x 10~^. We observe that if 
wc cut off the calculations for more lower values of survival probability it does not affect 
MFPT but the variances are slightly affected. Also if one approximates the sinusoidal signal 
better, by more than a six-step periodic signal one could obtain the positions of the minima 
of variances at exactly the same places as those with MFPT. 

From Fig. 9 we observe a slight deviation of the peak positions but we believe that if the 
sinusoidal signal is approximated by more than a six-step periodic signal the position of all 
the peaks will be the same. 
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FIGURES 



FIG. 1. Sinusoidal signal(dashed curve) and approximated six-step {p = 5) periodic signal(solid 
curve) for the full one cycle as a function of 6. 

FIG. 2. MFPT < d > as a function of u ; (a)for p=0 ;tlie usual telegraph signal (b) for p=l ; 
the two-step periodic signal, [A = 20,^o = A/2]. 

FIG. 3. MFPT < e{uj) > as a function of frequency w ; (a)A = 10, (b)A = 20, (c)A = 30, 

(d) A = 40, (e)A = 50, [p = 5,^0 = A/2]. 

FIG. 4. The variance as a function of co ; (a)A = 10, (b)A = 20, (c)A = 30, (d)A = 40, 

(e) A = 50, [p = 5,6 = A/2]. 

FIG. 5. MFPT < 6* > as a function of /? for length A = 30 ; (a) for resonant frequncy w* = 0.06 
(b) for off-resonant frequncy a; = 0.1 (c) for off-resonant frequncy u = 0.0, [p = 5]. 

FIG. 6. (a) FPTDF g{6) as a function of ^ for A = 30 before resonance for frequencies u = 
0.001, 0.003, 0.005, 0.01 (b) FPTDF g{e) as a function of 6* for A = 30 on and after resonance for 
frequencies lo = 0.06 (resonant), u = 0.07,0.09,0.1 respectively, [p = 5,^ = A/2]. 

FIG. 7. Resonant frequency u* as a function of length A. 

FIG. 8. MFPT at the resonant frequency < 0{uj*) > as a function of length A. 

FIG. 9. The dominant peaks of g{9)/uj* at resonant frequencies for different lengths 
(A = 20,25,30,35,40,50) as a function of lo*6. The lowermost curve is for A = 20 , and as 
length increases gradually upper curves are generated, [p = 5,^ = A/2]. 

FIG. 10. The universal functions a{x),lnU{x) as a function of x. 

FIG. 11. /a^x) as a function of x for A = 30 at resonance from actual data points (filled circle) 
[p = 5,6 = A/2] and fA{x) according to the formula (40) as a function of x for A = 30(solid 
curve). 



16 



0.012 



0.008 - 



0.004 - 



0.000 



0.012 - 




0.008 - 



0.004 - 



0.000 



T 1 r 

100 200 300 400 




0.012 



- 0.008 



- 0.004 



T 1 1 r 

50 100 150 200 250 



0.000 




T 1 1 r 

40 80 120 160 




0.020 



- 0.015 



- 0.010 



- 0.005 



0.000 



120 







1.2 
1.0 - 
0.8 - 
0.6 - 
0.4 - 
0.2 - 
0.0 



